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The problem of a particle confined in a box with moving walls is studied, focusing on the case of 
small perturbations which do not alter the shape of the boundary (‘pantography’). The presence 
of resonant transitions involving the natural transition frequencies of the system and the Fourier 
transform of the velocity of the walls of the box is brought to the light. The special case of a 
pantographic change of a circular box is analyzed in dept, also bringing to light the fact that the 
movement of the boundary cannot affect the angular momentum of the particle. 
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I. INTRODUCTION 

Resolution of Schrodinger equations with time-dependent Hamiltonian operators is challenging. In fact, exact 
resolutions are rare and limited to specific classes of problems PJ-Q - On the contrary, in most of the cases one 
can solve the dynamical problem only under special assumptions and with some approximations, as it happens for 
an Hamiltonian adiabatically changing [ 4 ] or in the presence of weak interactions which legitimate the use of a 
perturbative approach @-@1 • When the Hamiltonian time-dependence is periodical some special recipes based on the 
Floquet theory can be used 

In the panorama of systems with time-dependent Hamiltonians a special class is that of systems whose boundary 
conditions are time-dependent. Such kind of problems have been widely studied over the last decades in connection 
with the Casimir effect mm, but in the past it has also been considered in connection with quantum mechanical 
problems, from the Fermi quantum bouncer H to several works analyzing a “free” quantum particle in a box with 
moving walls jl5l - [l8j . The interest in such a class of problems is not only academic, and the relevance of moving 
boundaries has been discussed in cavity quantum electrodynamics El and in the physical scenario of trapped particles 
to propose new strategies for cooling atoms [20[. More recently, several other works appeared, dealing with specific box 
shapes [2ll . [22ll aimed at giving proper mathematical treatment of the Schrodinger problem in the presence of moving 
boundaries [ 23 , and exploring the raise of correlations between different particles confined in the same non-static [24j . 

In this paper, we analyze the dynamics of a quantum particle confined into a box whose walls are moving, and 
in particular we will focus on the possibility of exploiting moving boundaries to resonantly stimulate the system. 
Following the approach of our previous work [25| . we recast the original problem into the problem of a particle in a 
static box governed by a time-dependent Hamiltonian. The time-dependence of the boundary is then converted into a 
time-dependence of the relevant effective Hamiltonian acting in the fictitious static domain. Of course, if the walls move 
periodically, the effective time-dependent Hamiltonian will be periodic too, and with the same frequency. Through 
the exploitation of a time-dependent perturbative approach, we single out the presence of specific resonant transitions 
which is immediately traceable back to the analysis of the Fourier transform of the velocity of the boundary compared 
with the natural transition frequencies of the physical system. We will concentrate on two-dimensional problems 
but our results are immediately extendable to the three-dimensional case, and of course are essentially valid also for 
one-dimensional systems. 

The paper is organized as follows. In the next section we present the physical problem, the relevant mathematical 
formulation and the proper unitary transformation which removes the boundary motion and make it ‘replaced’ by a 
time-dependent term in the Hamiltonian. In the subsequent section we introduce the perturbative treatment singling 
out the presence of resonances and discussing the role of the Fourier transform of the boundary motion. In section HVl 
we specialize the previous results to the case of a two-dimensional circular box whose radius is oscillating at a single 
frequency. Finally, in section 0 we give some conclusive remarks. 


II. THE PHYSICAL SYSTEM AND HAMILTONIAN 
A. Two-dimensional box 

Let us consider a particle confined in a two-dimensional box whose contour, expressed in polar coordinates, is given 
by 


r = A(t)7(6»), d£[ 0,2tt], 


(1) 
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where A (t) is a smooth function of time, and 7 (0) is a smooth function of the angle 9 such that 7(27 t) = 7 ( 0 ) (which 
guarantees that the box is not open and the particle is confined). Because of the presence of A (t) the dimensions 
of the box change, but since the dependence of the radial coordinate of the box on 6 and t is factorized as X{t) r y{6) 
the shape of the box does not change. This is a kind of boundary modification that we call ‘pantographic’, and that 
corresponds to the fact that each point of the boundary moves only along the relevant ‘radial’ direction. 

The Hamiltonian describing the particle is simply given by the kinetic term: 

Tp fl 2 

H = 5- = -^, 

2 /i 2 /i 
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r dr ^ dr 2 r 2 d9 2 J ’ 

with the time-dependent Dirichlet boundary conditions: 

^(X(t) 7 (9),9) = 0. (3) 

More precisely, the Hamiltonian is given by an operator which coincides with the kinetic term in the domain 
corresponding to the box, and zero elsewhere. For the sake of simplicity we will avoid introduction of such a new 
symbol. 

According to our treatment in [23j and [25|, we make a unitary transformation which essentially maps the original 
domain delimited by the boundary given by Eq. © into another domain with the same shape but different diameter: 

r = 7 (9 ), 9 G [0, 2tt] . (4) 

We call the first domain (the time-dependent one) T>\ and this second one (the static one) V 1 . The relevant unitary 
operator acts as follows: 



( 2 a) 

( 2 b) 


(j){r, 9) = (U\ip)(r, 9) = Xi(>(rX, 9). 


( 5 ) 


It maps quantum states defined in T>\ into states defined in T> 1 . 

In the new picture, the operator which generates the dynamics of the particle is given by 


H e s 


U\HU\ + ihUxUt 


2/iA 2 


i h- 


1 + r 


dr 


( 6 ) 


where one has to take into account the fact that, since A depends on f, the unitary operator U\ depends on time as 
well. The Hamiltonian in Eq. © describes a particle with varying mass (because of the factor A “ 2 in the kinetic 
energy) in the presence of a time-dependent ‘potential’, i.e., the term \hX/X{\ + rd r ), which from now on we will call 
the ‘dilation potential’ or ‘dilation term’. 

It is remarkable to note that the dilation potential involves only the radial coordinate. This is a consequence of 
the pantographic nature of the box movement, which implies radial dilation and then only radial motion of the walls 
of the box. Of course, this simple fact does not guarantee conservation of the angular momentum. Indeed, generally 
speaking, the commutator between two operators depends on the domain in which the two operators are acting. 
Therefore, commutation between the dilation potential and the angular momentum will depend also on the shape of 
the box. 

It is also worth noting that in the more general case of changes with deformation (i.e., non pantographic) the 
effective Hamiltonian in the static domain turns out to be much more complicated than the one in Eq. © (the 
complete expression is reported in Ref. HBJ). 


B. One-dimensional and three-dimensional cases 


The one-dimensional and three-dimensional counterparts of our problem are treated in a very similar way. In the 
one-dimensional case, the domain is expressed as V\ = [—X(t)l/2,X(t)l/2] and will be mapped into = [—1/2,1/2], 
the boundary conditions are ip(±X(t)/2) = 0 and will be mapped into </>(±Z/2) = 0, and the generator of the time 
evolution will change as follows: 


H = - 


h 2 d 2 
2 /i dx 2 


H e s 


h 2 d 2 
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where the last differential operator can be put in the form rd r , with r = |x|. 

In the three-dimensional case, the domain is expressed as V\ = {(r, 0, p)\r < X(t)^(9, <p), 6 S [0, 2i r],tp £ [0, 7 r]} and 
will be mapped into D\ obtained for A = 1, the boundary conditions are ip(X(t)'y(9,ip),9,ip) = 0 V9,ip and will be 
mapped into cf>('y(9, ip), 0, ip) = 0, and the generator of the time evolution will change as follows: 


h 2 2 

H = -V 2 

2p 


-ffeff = - 


2pX 2 


X ( 3 d 
i h— —b r— 
A \ 2 dr 


( 8 ) 


This clarifies how the further results, which will be explicitly derived for the two-dimensional case, are essentially 
valid in ID and 3D. 


III. PANTOGRAPHIC PERTURBATIONS 

Once the problem of one particle in a 2D box with moving walls is transformed into the problem of a particle in a 
static box, we get the Hamiltonian in Eq. ([ 6 ]). 

Under the assumption that the dilation parameter A is a smooth and close-to-unity function, the quantity A/A is 
small and the dilation potential can be treated as a perturbation. This is in perfect agreement with considerations by 
Kato on the way to treat small modifications of the domain in which a particle can move [26| . Such an approach has 
already been exploited in Ref. [25| to treat changes of the domain with deformation. (It is appropriate to mention 
here that in [25| we assume knowledge of the dynamics in the pantographic case - the relevant Hamiltonian describing 
pantographic changes is considered as the unperturbed one- and then treat the terms coming from deformation as a 
perturbation. Nevertheless, knowledge of the dynamics in the pantographic case is quite limited, and exact dynamics 
is known only for the case of uniformly moving domain. In this paper, instead, we are attempting to extend our 
knowledge of the dynamics in the pantographic case.) 

Generally speaking, after making the unitary mapping to the domain in the new picture, that we address as 
the ‘Schrodinger picture’, we have a static domain and a time-dependent Hamiltonian (see Eq. (ftjjl l of the following 
form: 


H 3 (t) = X~ 2 Ho + XX-'V. (9) 

We express A as follows: 

A = 1 + ef(t), (10) 

where e is a dimensionless small parameter, f(t) is a smooth and bounded function, and the operators Hq and V are 
time-independent. We also introduce the series expansions of A -2 and A A” 1 : 

A 2 = ^4k e k = 1 — 2ef(t) + ..., ( 11 ) 

k 

A A - 1 = ^R k e k = e/(t) + ..., (12) 

k 

then getting the following Hamiltonian corrected to the first order in e: 

H s (t) = (l-2ef(t))H 0 + ef(t)V. (13) 

By performing the passage to the interaction picture through the unitary operator generated by the ‘unperturbed’ 
Hamiltonian X~ 2 Hq 7 Uo(t ) = e~ lh /o l A ( e > s )l we obtain the new generator of the time evolution, H\(t) = 

XX^U^VUoit). 

After performing the relevant spectral decompositions, 

Ho(t) = J2 E a \E a ){E a \ , (14) 

a 

and applying the standard perturbation treatment, we can explicitly write down the evolved state corrected up to the 
first order in the parameter e: 




e~Ti fo [l-2ef( g )]dsH 0 + e e --kH 0 t ^ 

cx.fi 



f(u)e^ E «- E ^ u du\E a )(E 0 \ 


\m) > 


(15) 
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where V a p = (E a \ V \Ep). 

On the basis of Eq. (1151) one can immediately argue that the movement of the walls induces quantum transitions, 
and two factors determine which transitions can occur: on the one hand the coupling through the dilation potential, 
which means that the matrix element V a p must be nonzero; on the other hand, in order to have a finite transition 
probability from the initial state to another at long time, the following condition must be satisfied: 


f(u) e^ Ea - E ^ u du + 0 , 


I E a — E f 3 1 


(16) 


Since the integral in Eq. [16] in the limit t — > oo essentially approaches the Fourier transform of f(t) for u; = —(E a — 
Ep)/h , one can assert that in order to have nonzero transitions between two states the Fourier transform of the radial 
velocity of the walls should be nonzero for the relevant frequency. On this basis, it is natural to talk about resonances 
and turns out to be important singling out the presence of resonant sinusoidal components in the perturbation (i.e., 
the dilation potential). Of course, as usual, non resonant components in the Fourier expansion of the perturbation 
can induce fluctuations determining small transitions at small time, but such transitions disappear after a sufficiently 
long time. 

As a relevant physical situation, one can consider the case where the A is a periodic function, and still smooth and 
close-to-unity. This means that f(t) is smooth and periodic, and so is f(t), easily leading to the fact that, at first 
order, only transitions associated to the frequency of / or its multiples are allowed. 


IV. THE ‘BREATHING’ CIRCLE 

In order to better fix the previous ideas we will consider a very special case, that we call the ‘breathing’ circle. 
In other words, we consider a particle moving inside a two-dimensional circular box whose radius is time-dependent: 
R(t) = Ro(l + esinwt), so that 7 (0) = Rq and /(f) = sinwt. After mapping the original problem into the problem of 
a particle confined in a circular box of radius Ro and expanding the relevant Hamiltonian with respect to e up to the 
first order, one obtains: 


H e g 


h 2 

-(1 — 2esinwt)V 2 few coswtifi 

2 /r 


dr 


(17) 


According to the analysis in Ref. [23], the eigenvalues and eigenfunctions of a ‘free’ particle (hence governed by 
Hq = — 7i 2 /(2^)V 2 ) in a 2D circular box are given by the following expressions: 


JR _ 

-C'ran r» 2 a mn ? 


2/irg 

Xmn = (2 tt) ^ X e , Tfl £ ^ , 

where J m (x) is the Bessel function of m-th order, a mn is the n-th zero of J mi and 


k 2 = 

mn 


2/rF' mn (u mn ) 


h 2 


fomn 
- 1/2 




JmiJ^mnX) dr 


(18) 

(19) 


( 20 ) 

( 21 ) 


We have already mentioned the fact that the operator 1 + rd r does not involve the angular variable but only the 
radial one. Moreover, because of the specific shape of the box, in this case the angular momentum of the system 
commutes with such an operator. This fact, considered also the commutation between kinetic energy and angular 
momentum, implies that the angular momentum is conserved through all the evolution. This result is indeed quite 
intuitive since the shape of the box is circular at any time and the ‘kicks’ that the particle can receive from the walls 
of the box are only directed along the radius of the circle, hence not producing any ‘torque’. On the contrary, because 
of such kicks, the energy of the particle can change, which corresponds to the possibility of having transitions between 
states with the same angular momentum but different energies. 

In connection with the perturbation treatment, we can say that the matrix elements of the dilation potential are 
vanishing when different values of the angular momentum quantum number are involved: V mnm / n r = (x mn | V |y m / ri /) = 
0 if m ^ m'. 
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FIG. 1: (Color online) Transition Probabilities to the states IX02) (blue bolded line), IX03) (red line) and |xo4) (black thin line) 
when the system is initially prepared in the ground state |xoi) and the frequency of the wall oscillations (<u) is equal to the 
frequency of the transition |xoi) —► |xo2). Here h = 1, fj, = 1, ro = 1 and e = 0.03. The time is expressed in units of 2n/co. 
Figures (a) and (b) show the probabilities in different time intervals. In the short-time scale (b) all the transitions are quite 
small. In the long-time scale (a) it is evident that the only transition probability which is non negligible is the one toward the 
state IX02), which is a clear manifestation of the resonance. 



FIG. 2: (Color online) The mean value of the distance of the particle from the centre of the box (r) when the system is prepared 
in its ground state |xoi) and the wall oscillation frequency uj is equal to the transition frequency from |xoi) to |xo2). Here 
h = 1, n = 1, ro = 1 and e = 0.03. The time is expressed in units of 2n/uj. 


In Fig. 03 we show the transition probabilities to some of the excited states when the system is prepared in the ground 
state and the oscillation of the walls resonates with the transition |xoi) —t |xo2). It is evident that the only significant 
transitions occur toward the state IX02), especially at relatively long time. The small and fast oscillation superimposed 
to the quadratic behavior of the population of IX02) is due to the fact that the rotating wave approximation has not 
been performed. 

In Fig. [^] we show how the mean value of the distance of the particle from the centre of the box ((r)) evolves when 
the particle is prepared in its ground states and the wall movement resonates with the transition |x'oi) —> |X02}- 


V. DISCUSSION 

The analysis developed in previous works, such as Ref. [23| and Ref. (2f|, show that when the problem of a system 
with moving boundaries is mapped into a problem with static boundaries the Hamiltonian inducing the evolution in 
the new picture becomes time-dependent. Even the very kinetic energy contains a time-dependent factor, so that the 
particle appears as a particle with changing mass. Moreover, a new term - the dilation term, which is time-dependent 
as well - is added. Such new term can induce transitions between the eigenstates of the unperturbed Hamiltonian, 
which in our case is simply given by the kinetic energy of a particle with a varying mass. Which transitions are 
induced depends on the way the boundary moves, and in particular from the Fourier transform of the velocity of the 
walls. As a very special case, we have considered the case where the walls oscillate at a precise frequency, so that 
transitions between eigenstates of the ‘free’ Hamiltonian can be induced when the boundary oscillations are properly 
‘tuned’. 

We emphasize the fact that our analysis demonstrates with an appropriate mathematical treatment the intuitive idea 
that a boundary oscillating at a given frequency can induce transitions, when such oscillations are tuned to a transition 
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frequency of the system. Then, in some sense, an oscillating boundary acts on the system like a suitable oscillating 
field. By this way, it is worth mentioning that, because of the choice of considering only pantographic changes of the 
domain, the dilation potential turns out to be a radial potential (1 +rd r ) and then, with an appropriate geometrical 
nature of the domain, it preserves angular momentum, as it happens in the case of a circular box. 

As a final remark, we point out that if the particle in the original time-dependent domain is subjected, not only 
to the potential describing the box contour, but also to another potential, of course such a term will be kept in the 
new picture with a suitable scaling (see Ref. [23}). In our case, for the sake of simplicity we have not considered 
such a situation. Nevertheless, in the presence of such additional potential the approach is essentially the same as 
before, just with some mathematical complications. In fact, in the static domain the particle would be described as 
a time varying particle subjected to a time-dependent potential (resulting as a time-dependent scaling of the original 
potential) and to the dilation term (resulting from the change of picture). Therefore, in the limit of small and smooth 
movements of the walls of the box, the last term will be treated as a perturbation to the time evolution induced by 
the first two terms. 
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